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Precise calculation of transition frequencies of hydrogen and deuterium 
based on a least-squares analysis 
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We combine a limited number of accurately measured transition frequencies in hydrogen and deuterium, 
recent quantum electrodynamics (QED) calculations, and, as an essential additional ingredient, a generalized 
least-squares analysis, to obtain precise and optimal predictions for hydrogen and deuterium transition frequen- 
cies. Some of the predicted transition frequencies have relative uncertainties more than an order of magnitude 
smaller than that of the g-factor of the electron, which was previously the most accurate prediction of QED. 

PACS numbers: 12.20.Ds, 31.30.Jv, 06.20.Jr, 31.15.-p 



In the past decade there have been significant advances in 
both the measurement and theory of transition frequencies in 
hydrogen and deuterium. With the advent of frequency-comb 
metrology, the progress in measurements has accelerated to 
the extent that the uncertainty in the experimental value of the 
1S-2S frequency in hydrogen has decreased by three orders 
of magnitude in about ten years 1 1]. Many other precision 
measurements of transitions in hydrogen and deuterium with 
principal quantum number n up to 12 have been made and are 
Usted in Table U together with the 1S-2S frequency flllll 
y, y, |6l tZl jsl y, lldl . In addition, precise measurements have 
been made involving states of hydrogen with n = 27 to 30 
lull . Further advances may be anticipated from a number of 
grovros currently working to improve measurement accuracy 

On the theoretical side, there has been progress in the accu- 
racy of quantum electrodynamics (QED) calculations which 
has led to a significant reduction of the uncertainty of both the 
one-photon 1 140111 El ESl and two-photon 1 2 1,0121 
124 l25l 12611 contributions. In addition to this progress, the ap- 
plication of the least-squares method described here makes it 
possible in some cases to predict transition frequencies with 
a relative uncertainty that is more than an order of magnitude 
smaller than the relative uncertainty in the Rydberg constant, 
which is 6.6 x 10~'^. 

It is of interest to have accurate calculations of hydrogen 
and deuterium transition frequencies for comparison to exist- 
ing and new experimental values, for both frequency standards 
applications, and as tests of QED. Also, with sufficiently ac- 
curate theory available, it would be possible to consider re- 
defining the SI second in terms of an assigned value of the 
Rydberg constant. Although this would entail about three or- 
ders of magnitude improvement in the theory, the recent rate 
of progress suggests that it may be within reach. 

In this letter, we give theoretical values for a number of 
transition frequencies that are meant to be best values consis- 
tent with cuiTently available experimental and theoretical in- 
formation. Such calculated values for the transitions in TableU 
follow from the results of the 2002 CODATA least-squares ad- 



justment of the fundamental constants 1271 l28l 12911 . and are 

listed in that table. Of course, the calculated value of the 
1S-2S transition in hydrogen should not be viewed as a theo- 
retical prediction, because the least-squares adjustment deter- 
mines values of the adjusted constants (variables or unknowns 
of the adjustment) so that the calculated frequency is essen- 
tially equal to the very precise measured value. The number 
in parentheses with each value is the standard uncertainty (es- 
timated standard deviation) in the last two figures. Hyperfine 
structure effects are not included in the quoted numbers. 

For hydrogen or deuterium transition frequencies not in- 
cluded in the 2002 adjustment, we proceed as follows. The 
energy level Ei of state i can be written as a function of 
the fundamental constants and an additional adjusted con- 
stant 5i which takes into account the uncertainty in the theory 
1 27, 30, 31]. For example, for the case in which i is a state of 
hydrogen, we have 



E, = Hi[R^,a,A,{e),A,{p),Rp\+5„ 



(1) 



where the constants that appear as arguments of the function 
Hi are listed in Tablellll Because the values of the constants in 
Eq. Q, including 5i, result from a least-squares adjustment, 
they are correlated, particularly those for i?oo and i?p, which 
have a correlation coefficient of 0.996. The uncertainty of 
the calculated value for the 1S-2S frequency in hydrogen is 
increased by a factor of about 500 if such correlations are ne- 
glected. The function Hi also depends implicitly on c and the 
Planck constant h. However, these constants are not displayed 
as variables, because c is a fixed number, and the frequencies 
{Ei — Ei')/h are essentially independent of h. Levels in deu- 
terium are given as similar functions with p replaced by d. 

The theory included in the function Hjjm Eq. Q is de- 
scribed in detail in Appendix A of Ref. [27^, which provides 
a review of the relevant calculations. Much of that informa- 
tion is in the form of equations that are valid for any state, 
with the exception of tables of data that only have entries for 
the levels included in the 2002 CODATA adjustment. En- 
larged versions of those tables with data for all states with 
n < 200 are available on the NIST Physics Laboratory Web 



TABLE I: Transition frequencies in iiydrogen uh and in deuterium i^u used in the 2002 CODATA least-squares adjustment of the values of 
the fundamental constants and the calculated values. Hyperfine effects are not included in these values. 



Experiment 


Frequency interval(s) 




Reported value 


Calculated value 












v/kRz 


///kHz 


Niering et al. ^ 


!^h(1Si/2 - 


-2S1/2) 






2 466 061413187.103(46) 


2 466 061413187.103(46) 


Weitz efaZ. [2] 


I^h(2Si/2 - 


-4S1/2)- 


i!^H(lSi/2 - 


-2S1/2) 


4 797 338(10) 


4 797 331.8(2.0) 




I^h(2Si/2 - 


-4D5/2)- 


3fH(lSi/2 


-2S1/2) 


6 490144(24) 


6 490129.9(1.7) 




J^d(2Si/2^ 


-4S1/2)- 


ii^D(lSl/2 - 


- 2S1/2) 


4 801693(20) 


4801710.2(2.0) 




I^d(2Si/2^ 


" 4D5/2) - 


- i!^D(lSl/2 


"2Si/2) 


6 494 841(41) 


6494 831.5(1.7) 


Huber et al. [^J 


J'd(1Si/2 - 


-2S1/2)- 


Wl(lSl/2 - 


2S1/2) 


670 994 334.64(15) 


670 994 334.64(15) 


de Beauvoir et al. [4] 


I^h(2Si/2 - 


-8S1/2) 






770 649 350 012.0(8.6) 


770 649 350 016.1(2.8) 




J^h(2Si/2 - 


- 8D3/2) 






770 649 504 450.0(8.3) 


770 649 504 449.1(2.8) 




I^h(2Si/2 ' 


- 8D5/2) 






770 649 561584.2(6.4) 


770 649 561578.2(2.8) 




I^d(2Si/2 - 


- 8S1/2) 






770 859 041245.7(6.9) 


770 859 041242.6(2.8) 




1^d(2Si/2 - 


- 8D3/2) 






770 859195 701.8(6.3) 


770 859195 700.3(2.8) 




1^d(2Si/2 - 


- 8D5/2) 






770 859 252 849.5(5.9) 


770 859 252 845.1(2.8) 


Schwob et al. [5J 


i/h(2Si/2 - 


- 12D3/2) 






799191710 472.7(9.4) 


799191710 481.9(3.0) 




I^h(2Si/2 - 


- 12D5/2) 






799191727403.7(7.0) 


799191727409.1(3.0) 




i^D(2Si/2 - 


- 12D3/2) 






799 409168 038.0(8.6) 


799 409168 041.7(3.0) 




1^0(281/2^ 


- 12D5/2) 






799 409184 966.8(6.8) 


799 409184 973.4(3.0) 


Bourzeix et al. \6^ 


I^h(2Si/2 - 


-6S1/2)- 


|!^h(1Si/2 ^ 


- 3S1/2) 


4197 604(21) 


4197 600.3(2.2) 




I^h(2Si/2 - 


- 6D5/2) - 


i!^H(lSi/2 


— 3S1/2) 


4 699 099(10) 


4 699105.4(2.2) 


Berkeland et al. [7] 


J^h(2Si/2 - 


-4P1/2)- 




-2S1/2) 


4 664 269(15) 


4 664 254.3(1.7) 




I^h(2Si/2 - 


-4P3/2)- 


-2S1/2) 


6 035 373(10) 


6035 384.1(1.7) 


Hagley and Pipkin [M 


I^h(2Si/2 - 


-2P3/2) 






9 911200(12) 


9911197.6(2.4) 


Lundeen and Pipkin [g] 


I/h(2Pi/2 


-2S1/2) 






1057 845.0(9.0) 


1057 844.0(2.4) 


Newton etal. [10] 


i^H(2Pi/2 


-2S1/2) 






1057 862(20) 


1057 844.0(2.4) 



site at physics.nist.gov/hdel. Estimates of the theoretical un- 
certainties of the function Hi, represented by the constant Si 
in Eq. Q, are also given in Appendix A of Ref. 127). The a 
priori estimated value of 6i is Si (th) = 0, because the the- 
oretical expression for the levels includes all known contri- 
butions. However, the estimated uncertainty w[(5i(th)] is not 
zero, and there are significant covariances between the vari- 
ous (5s that take into account the expected patterns in the un- 
certainties. For example, for S states there are components 
of uncertainty with the functional form C jr? , where C is a 
common unknown constant, and there are components of un- 
certainty common to hydrogen and deuterium levels with the 
same quantum numbers. The theoretical uncertainties and co- 
variances are included in the least-squares adjustment as input 

TABLE II: The CODATA 2002 values of the constants used in the 
evaluation of the spectrum of hydrogen and deuterium. 



Constant 



Value 



Speed of light 
Rydberg constant 
Fine-structure constant 
Electron relative atomic mass 
Proton relative atomic mass 
Deuteron relative atomic mass 
Proton rms charge radius 
Deuteron rms charge radius 



c = 299 792 458ms~^ 

R^^ 10 973 731.568 525(73) m'^ 

0=1/137.035 99911(46) 

ylr(e) = 5.485 799 0945(24)xl0-'' 

^r(p) = 1.007276 466 88(13) 

Ar (d) = 2.013 553 212 70(35) 

7?p = 0.8750(68) fm 

fld = 2. 1394(28) fm 



data for the adjusted variables bi. 

A generalized least-squares adjustment is formulated here 
along the lines described in Refs. ll29ll and 12811 . Symbols that 
refer to data used in the 2002 CODATA adjustment of the con- 
stants are also defined in Ref. 1.2 8.1 . New energy levels Ei to 
be determined are added to the adjustment, along with the cor- 
responding theoretical expressions of the form in Eq. 0, and 
for each added level not among those in Table HJ a new ad- 
justed variable (5; is added. The updated column vector of 
input data Qu, matrix of their covariances Vu, and column 
vector of variables Z^ are written in block form as 



Qu 



Q \ 




V T 


Qs 


; K = 


T^ S 


Qe) 




\ Ve 



z„ 



Zs 



, (2) 



where Q, V, and Z are the corresponding sets of quantities 
used in the 2002 least-squares adjustment, Qs is the set of 
theoretical data 6i{th) ~ for the new variables Si, Zs is the 
new set of adjusted variables Si, and Qe is input data for the 
new energy levels Ei. In K, where Kifc = cov{Qui,Quk), 
S and T are the sets of theoretical covariances involving the 
new Ss, and Ve is the set of covariances of the new levels 
El. Since the input data for the new levels are unknown, we 
simply assume that the uncertainties are very large and that 
there are no correlations among them or with the rest of the 
input data. This yields the blocks of zeros in Vu and results in 



Ve being diagonal. 

The input data and adjusted variables are related by the set 
of observational equations given by 



Qu — Fu{Zu); 





(3) 



where the dot over the equal sign indicates that the equa- 
tion represents the ideal relations between the input data and 
the adjusted constants which are not simultaneously satisfied, 
since the set of equations is overdetermined. In Eq. (|3}, F is 
the set of functions in the observational equations of the 2002 
adjustment, and E is the set of expressions for the new energy 
levels of the form in Eq. Q. The observational equations are 
linearized by writing the Taylor series 



)u = Fu(Z(0)) + A„(Z„-Z(")) + ..., 



where A,, is the matrix of derivatives 



aK„(zi' 



(0)n 



^UlJ 



dZ, 



(0) 



A, 




(4) 



(5) 



and neglecting higher-order terms. In Eq. (|5|l, A is the ma- 
trix of derivatives from the 2002 adjustment, / is the identity 
matrix, and B and C are derivatives of the new energy levels 
with respect to the old and new variables, respectively. The 
truncated expression in Eq. (|4} corresponds to 



where Yu - Qu - Fu{Z^°^) and Xu 
The update adjustment starts with 



(6) 



Zu 



Z^\ 



Z(") 



(7) 



where Z is the final vector of constants from the 2002 adjust- 
ment and 



i^, 




Qs 



F{Z) 
7(0) 



(8) 



Qe - E(Z^°^^ 



where Y is the final value of Y from the 2002 adjustment and 
Ys = 0. The solution X^ to Eq. (|6|l, the value of X^ that 

minimizes (Y^ - A^X^yV^^^iYu - ^u-'^u), is 

lu = G^A'^V-'Y^; Gu = (A'^V-^A^yK (9) 

The covariance matrix of the solution Xu is Gu, and its cal- 
culation is the key to the update. The Schur-Banachiewicz 
inverse formula |32, 33] applied to the upper-left four blocks 
of the matrix Vu in Eq. (|2j gives 

''V-^+V-^TRT^V-^ -V-^TR 
V-^ = I -RT^V-^ R 1, (10) 











vz 



where R = {S — T^V ^T) ^. For increasing uncertain- 
ties of the unknown input data for the new levels Ei, we have 



V„ 



0, and we work in this limit. A direct calculation 



from Eqs. (|5}, (|9}, and ^, with V^^ = 0, yields 



G,, — 



G-^ + G-^URU^G-^ 
-RU^G"^ 



-G-^UR 
R 



(11) 



where G = {A^V^^Ay^ is the matrix from the 2002 ad- 
justment and U = GA^V^^T. Evidently, Eq. (E) is the 
Schur-Banachiewicz inverse expression for 



Gu — 



G U 

U^ P 



(12) 



provided R = {P - U^G-'^U)-'^, that is, if 

P = S -T^V^^T + U'^G-^U ^ S + DT, (13) 

where D = T^y-^ {AGA'^V-^ - l). This result for Gu in 
terms of G means that the exact result of the enlarged least- 
squares adjustment can be obtained from results of the 2002 
least-squares adjustment with a relatively simple calculation. 
That is, the matrix inversions needed for the enlarged adjust- 
ment have effectively been carried out exactly, with the results 
explicitly expressed in terms of the matrices and vectors of the 
2002 adjustment. In particular. 



GuAjV-' 



GA'^V-^ 
D I 



so that 

^u = 

or for the adjusted constants 



GA'^V-^lQ-FiZ)] \ 
D[Q~F{Z)] J 




DY 



Zu = 4"^ + ^u 



z 

DY 



(14) 



(15) 



(16) 



with covariance matrix cov(Zu) — Gu- More importantly, 
Eqs. J12t and ( I16t show that both the values and uncertain- 
ties of the new levels being calculated are influenced by their 
covariances with the levels in the 2002 least-squares adjust- 
ment, while the values and uncertainties of the variables from 
that adjustment are not changed at all. Also, since the only 
adjusted variables that change in the update appear linearly in 
Eq. (|3j, no iteration of the update is needed to reach the final 
result. 

The energy levels and their covariances are thus given by 



Qi 



E Z„ 



(17) 



cov(Qb) = BGB^ + CU^B^ +BUC^ + CPC^ 



where the latter result is the lower-right block of the relation 
cov{Qu) = ^uGu^u . The result from Eq. ( I17t for a transi- 
tion frequency i/im and its standard uncertainty u{i'im) for the 



TABLE III: Calculated transition frequencies in hydrogen and deu- 
terium from the IS state to the 3S and 3D excited states. 



Excited 
state 



Hydrogen 
i^H/kHz 



Deuterium 



3Si/2 2 922 743 278 671.6(1.4) 
3D3/2 2 922 746 208 551.40(70) 
3D5/2 2 922 747 291888.61(70) 



2 923 538 534 391.8(1.4) 
2 923 541464 741.75(72) 
2 923 542 548 374.66(72) 



TABLE IV: Examples of calculated transition frequencies in hydro- 
gen and deuterium from the 2S state to various S and D excited states. 



Excited 


Hydrogen 


Deuterium 


state 


i^H/kHz 


j/D/kHz 


3Si/2 


456 681865 484.5(1.4) 


456 806126 870.1(1.4) 


3D3/2 


456 684 795 364.30(69) 


456 809 057 220.01(69) 


3D5/2 


456 685 878 701.51(69) 


456 810140 852.91(69) 


4Si/2 


616 520150 628.5(2.0) 


616 687 903 590.7(2.0) 


4D3/2 


616 521386 393.3(1.7) 


616 689139 553.8(1.7) 


4D5/2 


616 521843 426.7(1.7) 


616 689 596 711.9(1.7) 



transition I ^ mis given by 

hl^lm = QeI — QEm (18) 

1 

hu{vtm) = u'^{Qei) -'2cOv{QEhQEm) +u'^{QE7n) 

where u^{Qei) = cov{Qei, Qei), i ^ l,m. 

Examples of calculated transition frequencies in hydrogen 
and deuterium based on this update, starting from the results 
of the 2002 least-squares adjustment, are given in Tables Hill 
and II VI Data from that adjustment needed for such a calcu- 
lation are available on the Web at physics.nist.gov/constants. 
The frequencies in Tables Illlland llVl all have relative uncer- 
tainties that are smaller than the relative uncertainty of the 
Rydberg constant, mainly as a result of the correlations be- 
tween Roo, Rp, and i?d- In some cases, these values are nearly 
five orders of magnitude more accurate than the correspond- 
ing best previous values |34]. A database that gives the fre- 
quency of any transition between levels with n < 200 based 
on the calculations described here is maintained on the Web 
at physics.nist.gov/hdel. 

Helpful conversations with G. W. Stewart are acknowl- 
edged by one of the authors (PJM). 
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